MicroRNAs modulate immunological and inflammatory responses in Holstein cattle naturally infected with Mycobacterium avium subsp. paratuberculosis.

MicroRNAs (miRNAs) regulate the post-transcriptional expression of genes by binding to their target mRNAs. In this study, whole miRNA sequencing was used to compare the expression of miRNAs in ileocecal valve (ICV) and peripheral blood (PB) samples of cows with focal or diffuse paratuberculosis (PTB)-associated lesions in gut tissues versus (vs) control cows without lesions. Among the eight miRNAs differentially expressed in the PB samples from cows with diffuse lesions vs controls, three (miR-19a, miR-144, miR32) were also down-regulated in cows with diffuse vs focal lesions. In the ICV samples, we identified a total of 4, 5, and 18 miRNAs differentially expressed in cows with focal lesions vs controls, diffuse lesions vs controls, and diffuse vs focal lesions, respectively. The differential expression of five microRNAs (miR-19a, miR-144, miR-2425-3p, miR-139, miR-101) was confirmed by RT-qPCR. Next, mRNA target prediction was performed for each differentially expressed miRNA. A functional analysis using the predicted gene targets revealed a significant enrichment of the RNA polymerase and MAPK signaling pathways in the comparison of cows with focal vs no lesions and with diffuse vs focal lesions, respectively. The identified miRNAs could be used for the development of novel diagnostic and therapeutical tools for PTB control.

Multi-omics Characterization of Response to PD-1 Inhibitors in Advanced Melanoma.

Immunotherapy improves the survival of patients with advanced melanoma, 40% of whom become long-term responders. However, not all patients respond to immunotherapy. Further knowledge of the processes involved in the response and resistance to immunotherapy is still needed. In this study, clinical paraffin samples from fifty-two advanced melanoma patients treated with anti-PD-1 inhibitors were assessed via high-throughput proteomics and RNA-seq. The obtained proteomics and transcriptomics data were analyzed using multi-omics network analyses based on probabilistic graphical models to identify those biological processes involved in the response to immunotherapy. Additionally, proteins related to overall survival were studied. The activity of the node formed by the proteins involved in protein processing in the endoplasmic reticulum and antigen presentation machinery was higher in responders compared to non-responders; the activity of the immune and inflammatory response node was also higher in those with complete or partial responses. A predictor for overall survival based on two proteins (AMBP and PDSM5) was defined. In summary, the response to anti-PD-1 therapy in advanced melanoma is related to protein processing in the endoplasmic reticulum, and also to genes involved in the immune and inflammatory responses. Finally, a two-protein predictor can define survival in advanced disease. The molecular characterization of the mechanisms involved in the response and resistance to immunotherapy in melanoma leads the way to establishing therapeutic alternatives for patients who will not respond to this treatment.

Comparison of Extracellular Vesicle Isolation Methods for miRNA Sequencing.

MicroRNAs (miRNAs) encapsulated in extracellular vesicles (EVs) are potential diagnostic and prognostic biomarkers. However, discrepancies in miRNA patterns and their validation are still frequent due to differences in sample origin, EV isolation, and miRNA sequencing methods. The aim of the present study is to find a reliable EV isolation method for miRNA sequencing, adequate for clinical application. To this aim, two comparative studies were performed in parallel with the same human plasma sample: (i) isolation and characterization of EVs obtained using three procedures: size exclusion chromatography (SEC), iodixanol gradient (GRAD), and its combination (SEC+GRAD) and (ii) evaluation of the yield of miRNA sequences obtained using NextSeq 500 (Illumina) and three miRNA library preparation protocols: NEBNext, NEXTFlex, and SMARTer smRNA-seq. The conclusion of comparison (i) is that recovery of the largest amount of EVs and reproducibility were attained with SEC, but GRAD and SEC+GRAD yielded purer EV preparations. The conclusion of (ii) is that the NEBNext library showed the highest reproducibility in the number of miRNAs recovered and the highest diversity of miRNAs. These results render the combination of GRAD EV isolation and NEBNext library preparation for miRNA retrieval as adequate for clinical applications using plasma samples.

Client Applications and Server-Side Docker for Management of RNASeq and/or VariantSeq Workflows and Pipelines of the GPRO Suite.

The GPRO suite is an in-progress bioinformatic project for -omics data analysis. As part of the continued growth of this project, we introduce a client- and server-side solution for comparative transcriptomics and analysis of variants. The client-side consists of two Java applications called "RNASeq" and "VariantSeq" to manage pipelines and workflows based on the most common command line interface tools for RNA-seq and Variant-seq analysis, respectively. As such, "RNASeq" and "VariantSeq" are coupled with a Linux server infrastructure (named GPRO Server-Side) that hosts all dependencies of each application (scripts, databases, and command line interface software). Implementation of the Server-Side requires a Linux operating system, PHP, SQL, Python, bash scripting, and third-party software. The GPRO Server-Side can be installed, via a Docker container, in the user's PC under any operating system or on remote servers, as a cloud solution. "RNASeq" and "VariantSeq" are both available as desktop (RCP compilation) and web (RAP compilation) applications. Each application has two execution modes: a step-by-step mode enables each step of the workflow to be executed independently, and a pipeline mode allows all steps to be run sequentially. "RNASeq" and "VariantSeq" also feature an experimental, online support system called GENIE that consists of a virtual (chatbot) assistant and a pipeline jobs panel coupled with an expert system. The chatbot can troubleshoot issues with the usage of each tool, the pipeline jobs panel provides information about the status of each computational job executed in the GPRO Server-Side, while the expert system provides the user with a potential recommendation to identify or fix failed analyses. Our solution is a ready-to-use topic specific platform that combines the user-friendliness, robustness, and security of desktop software, with the efficiency of cloud/web applications to manage pipelines and workflows based on command line interface software.

Formal Meta-Analysis of Hypoxic Gene Expression Profiles Reveals a Universal Gene Signature.

Integrating transcriptional profiles results in identifying gene expression signatures that are more robust than those obtained for individual datasets. However, a direct comparison of datasets derived from heterogeneous experimental conditions is problematic, hence their integration requires applying of specific meta-analysis techniques. The transcriptional response to hypoxia has been the focus of intense research due to its central role in tissue homeostasis and prevalent diseases. Accordingly, many studies have determined the gene expression profile of hypoxic cells. Yet, despite this wealth of information, little effort has been made to integrate these datasets to produce a robust hypoxic signature. We applied a formal meta-analysis procedure to datasets comprising 430 RNA-seq samples from 43 individual studies including 34 different cell types, to derive a pooled estimate of the effect of hypoxia on gene expression in human cell lines grown ingin vitro. This approach revealed that a large proportion of the transcriptome is significantly regulated by hypoxia (8556 out of 20,888 genes identified across studies). However, only a small fraction of the differentially expressed genes (1265 genes, 15%) show an effect size that, according to comparisons to gene pathways known to be regulated by hypoxia, is likely to be biologically relevant. By focusing on genes ubiquitously expressed, we identified a signature of 291 genes robustly and consistently regulated by hypoxia. Overall, we have developed a robust gene signature that characterizes the transcriptomic response of human cell lines exposed to hypoxia in vitro by applying a formal meta-analysis to gene expression profiles.

SARS-CoV-2 Mutant Spectra at Different Depth Levels Reveal an Overwhelming Abundance of Low Frequency Mutations.

Populations of RNA viruses are composed of complex and dynamic mixtures of variant genomes that are termed mutant spectra or mutant clouds. This applies also to SARS-CoV-2, and mutations that are detected at low frequency in an infected individual can be dominant (represented in the consensus sequence) in subsequent variants of interest or variants of concern. Here we briefly review the main conclusions of our work on mutant spectrum characterization of hepatitis C virus (HCV) and SARS-CoV-2 at the nucleotide and amino acid levels and address the following two new questions derived from previous results: (i) how is the SARS-CoV-2 mutant and deletion spectrum composition in diagnostic samples, when examined at progressively lower cut-off mutant frequency values in ultra-deep sequencing; (ii) how the frequency distribution of minority amino acid substitutions in SARS-CoV-2 compares with that of HCV sampled also from infected patients. The main conclusions are the following: (i) the number of different mutations found at low frequency in SARS-CoV-2 mutant spectra increases dramatically (50- to 100-fold) as the cut-off frequency for mutation detection is lowered from 0.5% to 0.1%, and (ii) that, contrary to HCV, SARS-CoV-2 mutant spectra exhibit a deficit of intermediate frequency amino acid substitutions. The possible origin and implications of mutant spectrum differences among RNA viruses are discussed.

Hypoxia classifier for transcriptome datasets.

Molecular gene signatures are useful tools to characterize the physiological state of cell populations, but most have developed under a narrow range of conditions and cell types and are often restricted to a set of gene identities. Focusing on the transcriptional response to hypoxia, we aimed to generate widely applicable classifiers sourced from the results of a meta-analysis of 69 differential expression datasets which included 425 individual RNA-seq experiments from 33 different human cell types exposed to different degrees of hypoxia (0.1-5%[Formula: see text]) for 2-48 h. The resulting decision trees include both gene identities and quantitative boundaries, allowing for easy classification of individual samples without control or normoxic reference. Each tree is composed of 3-5 genes mostly drawn from a small set of just 8 genes (EGLN1, MIR210HG, NDRG1, ANKRD37, TCAF2, PFKFB3, BHLHE40, and MAFF). In spite of their simplicity, these classifiers achieve over 95% accuracy in cross validation and over 80% accuracy when applied to additional challenging datasets. Our results indicate that the classifiers are able to identify hypoxic tumor samples from bulk RNAseq and hypoxic regions within tumor from spatially resolved transcriptomics datasets. Moreover, application of the classifiers to histological sections from normal tissues suggest the presence of a hypoxic gene expression pattern in the kidney cortex not observed in other normoxic organs. Finally, tree classifiers described herein outperform traditional hypoxic gene signatures when compared against a wide range of datasets. This work describes a set of hypoxic gene signatures, structured as simple decision tress, that identify hypoxic samples and regions with high accuracy and can be applied to a broad variety of gene expression datasets and formats.

miRNA Profile Based on ART Delay in Vertically Infected HIV-1 Youths Is Associated With Inflammatory Biomarkers and Activation and Maturation Immune Levels.

Early antiretroviral treatment (ART) in vertically acquired HIV-1-infection is associated with a rapid viral suppression, small HIV-1 reservoir, reduced morbimortality and preserved immune functions. We investigated the miRNA profile from vertically acquired HIV-1-infected young adults based on ART initiation delay and its association with the immune system activation. Using a microRNA panel and multiparametric flow cytometry, miRNome profile obtained from peripheral blood mononuclear cells and its association with adaptive and innate immune components were studied on vertically HIV-1-infected young adults who started ART early (EARLY, 0-53 weeks after birth) and later (LATE, 120-300 weeks). miR-1248 and miR-155-5p, were significantly upregulated in EARLY group compared with LATE group, while miR-501-3p, miR-548d-5p, miR-18a-3p and miR-296-5p were significantly downregulated in EARLY treated group of patients. Strong correlations were obtained between miRNAs levels and soluble biochemical biomarkers and immunological parameters including CD4 T-cell count and maturation by CD69 expression on CD4 T-cells and activation by HLA-DR on CD16high NK cell subsets for miR-1248 and miR-155-5p. In this preliminary study, a distinct miRNA signature discriminates early treated HIV-1-infected young adults. The role of those miRNAs target genes in the modulation of HIV-1 replication and latency may reveal new host signaling pathways that could be manipulated in antiviral strategies. Correlations between miRNAs levels and inflammatory and immunological markers highlight those miRNAs as potential biomarkers for immune inflammation and activation in HIV-1-infected young adults who initiated a late ART.

Vaccine breakthrough infections with SARS-CoV-2 Alpha mirror mutations in Delta Plus, Iota, and Omicron.

Replication of SARS-CoV-2 in the human population is defined by distributions of mutants that are present at different frequencies within the infected host and can be detected by ultra-deep sequencing techniques. In this study, we examined the SARS-CoV-2 mutant spectra of amplicons from the spike-coding (S-coding) region of 5 nasopharyngeal isolates derived from patients with vaccine breakthrough. Interestingly, all patients became infected with the Alpha variant, but amino acid substitutions that correspond to the Delta Plus, Iota, and Omicron variants were present in the mutant spectra of the resident virus. Deep sequencing analysis of SARS-CoV-2 from patients with vaccine breakthrough revealed a rich reservoir of mutant types and may also identify tolerated substitutions that can be represented in epidemiologically dominant variants.

Description of the genetic variants identified in a cohort of patients diagnosed with localized anal squamous cell carcinoma and treated with panitumumab.

Squamous cell carcinoma is the most frequent histologic type of anal carcinoma. The standard of care since the 1970s has been a combination of 5-fluorouracil, mitomycin C, and radiotherapy. This treatment is very effective in T1/T2 tumors (achieving complete regression in 80-90% of tumors). However, in T3/T4 tumors, the 3-year relapse free survival rate is only 50%. The VITAL trial aimed to assess the efficacy and safety of panitumumab in combination with this standard treatment. In this study, 27 paraffin-embedded samples from the VITAL trial and 18 samples from patients from daily clinical practice were analyzed by whole-exome sequencing and the influence of the presence of genetic variants in the response to panitumumab was studied. Having a moderate- or high-impact genetic variant in PIK3CA seemed to be related to the response to panitumumab. Furthermore, copy number variants in FGFR3, GRB2 and JAK1 were also related to the response to panitumumab. These genetic alterations have also been studied in the cohort of patients from daily clinical practice (not treated with panitumumab) and they did not have a predictive value. Therefore, in this study, a collection of genetic alterations related to the response with panitumumab was described. These results could be useful for patient stratification in new anti-EGFR clinical trials.

Immune expression profile identification in a group of proliferative verrucous leukoplakia patients: a pre-cancer niche for oral squamous cell carcinoma development.

OBJECTIVES: To explore the pathophysiology of proliferative verrucous leukoplakia, a rare oral disorder that exhibits high rates of recurrence and malignant transformation, through a RNAseq case-control study. MATERIAL AND METHODS: We obtained oral biopsies from 10 patients with verrucous leukoplakia lesions and from the mucosa of 5 healthy individuals for sequencing using RNAseq technology. Using bioinformatic methods, we investigated gene expression and enrichment differences between patients both with and without the disorder. We applied network biology methods to investigate functional relations among those genes that were differentially deregulated. RESULTS: We detected 140 differentially expressed genes with distinct roles in immune surveillance, tissue and organ morphogenesis, development, and organization. Of these 140 genes, 111 have been previously described as cancer expression biomarkers, being oral squamous cell carcinoma the most represented type of cancer among them. Of these 140 genes, 26 were prioritized for further investigation as biomarkers using larger sample sizes. CONCLUSIONS: The gene expression patterns of healthy and unhealthy patients differed in 140 genes whose deregulation has a functional impact on normal functioning of the immune system. This immune expression profile provides a plausible hypothesis to explain the transformation to oral squamous cell carcinoma observed in 6 of the 10 assayed cases. CLINICAL RELEVANCE: By determining the molecular bases of the proliferative verrucous leukoplakia disorder and identifying early biomarkers of malignancy, this can allow us to develop new treatment strategies.

Validation and clinical application of a targeted next-generation sequencing gene panel for solid and hematologic malignancies.

BACKGROUND: Next-generation sequencing (NGS) is a high-throughput technology that has become widely integrated in molecular diagnostics laboratories. Among the large diversity of NGS-based panels, the Trusight Tumor 26 (TsT26) enables the detection of low-frequency variants across 26 genes using the MiSeq platform. METHODS: We describe the inter-laboratory validation and subsequent clinical application of the panel in 399 patients presenting a range of tumor types, including gastrointestinal (GI, 29%), hematologic (18%), lung (13%), gynecological and breast (8% each), among others. RESULTS: The panel is highly accurate with a test sensitivity of 92%, and demonstrated high specificity and positive predictive values (95% and 96%, respectively). Sequencing testing was successful in two-thirds of patients, while the remaining third failed due to unsuccessful quality-control filtering. Most detected variants were observed in the TP53 (28%), KRAS (16%), APC (10%) and PIK3CA (8%) genes. Overall, 372 variants were identified, primarily distributed as missense (81%), stop gain (9%) and frameshift (7%) altered sequences and mostly reported as pathogenic (78%) and variants of uncertain significance (19%). Only 14% of patients received targeted treatment based on the variant determined by the panel. The variants most frequently observed in GI and lung tumors were: KRAS c.35G > A (p.G12D), c.35G > T (p.G12V) and c.34G > T (p.G12C). CONCLUSIONS: Prior panel validation allowed its use in the laboratory daily practice by providing several relevant and potentially targetable variants across multiple tumors. However, this study is limited by high sample inadequacy rate, raising doubts as to continuity in the clinical setting.

Altered DNA methylation in human placenta after (suspected) preterm labor.

Aim: The aim of this study was to determine if alterations in DNA methylation in the human placenta would support suspected preterm labor as a pathologic insult associated with diminished placental health. Methods: We evaluated placental DNA methylation at seven loci differentially methylated in placental pathologies using targeted bisulfite sequencing, in placentas associated with preterm labor (term birth after suspected preterm labor [n = 15] and preterm birth [n = 15]), and controls (n = 15). Results: DNA methylation levels at the NCAM1 and PLAGL1 loci in placentas associated with preterm labor did differ significantly (p < 0.05) from controls. Discussion: Specific alterations in methylation patterns indicative of an unfavourable placental environment are associated with preterm labor per se and not restricted to preterm birth.

Increased miR-7641 Levels in Peritoneal Hyalinizing Vasculopathy in Long-Term Peritoneal Dialysis Patients.

Peritoneal hyalinizing vasculopathy (PHV) represents the cornerstone of long-term peritoneal dialysis (PD), and especially characterizes patients associated with encapsulating peritoneal sclerosis. However, the mechanisms of PHV development remain unknown. A cross sectional study was performed in 100 non-selected peritoneal biopsies of PD patients. Clinical data were collected and lesions were evaluated by immunohistochemistry. In selected biopsies a microRNA (miRNA)-sequencing analysis was performed. Only fifteen patients (15%) showed PHV at different degrees. PHV prevalence was significantly lower among patients using PD fluids containing low glucose degradation products (GDP) (5.9% vs. 24.5%), angiotensin converting enzyme inhibitors (ACEIs) (7.5% vs. 23.4%), statins (6.5% vs. 22.6%) or presenting residual renal function, suggesting the existence of several PHV protective factors. Peritoneal biopsies from PHV samples showed loss of endothelial markers and induction of mesenchymal proteins, associated with collagen IV accumulation and wide reduplication of the basement membrane. Moreover, co-expression of endothelial and mesenchymal markers, as well as TGF-ß1/Smad3 signaling activation were found in PHV biopsies. These findings suggest that an endothelial-to-mesenchymal transition (EndMT) process was taking place. Additionally, significantly higher levels of miR-7641 were observed in severe PHV compared to non-PHV peritoneal biopsies. Peritoneal damage by GDPs induce miRNA deregulation and an EndMT process in submesothelial vessels, which could contribute to collagen IV accumulation and PHV.

Comprehensive Characterization of the Mutational Landscape in Localized Anal Squamous Cell Carcinoma.

Anal squamous cell carcinoma (ASCC) is a rare neoplasm. Chemoradiotherapy is the standard of care, with no therapeutic advances achieved over the past three decades. Thus, a deeper molecular characterization of this disease is still necessary. We analyzed 46 paraffin-embedded tumor samples from patients diagnosed with primary ASCC by exome sequencing. A bioinformatics approach focused in the identification of high-impact genetic variants, which may act as drivers of oncogenesis, was performed. The relation between genetics variants and prognosis was also studied. The list of high-impact genetic variants was unique for each patient. However, the pathways in which these genes are involved are well-known hallmarks of cancer, such as angiogenesis or immune pathways. Additionally, we determined that genetic variants in BRCA2, ZNF750, FAM208B, ZNF599, and ZC3H13 genes are related with poor disease-free survival in ASCC. This may help to stratify the patient's prognosis and open new avenues for potential therapeutic intervention. In conclusion, sequencing of ASCC clinical samples appears an encouraging tool for the molecular portrait of this disease.

Genetic Profile and Functional Proteomics of Anal Squamous Cell Carcinoma: Proposal for a Molecular Classification.

Anal squamous cell carcinoma is a rare tumor. Chemo-radiotherapy yields a 50% 3-year relapse-free survival rate in advanced anal cancer, so improved predictive markers and therapeutic options are needed. High-throughput proteomics and whole-exome sequencing were performed in 46 paraffin samples from anal squamous cell carcinoma patients. Hierarchical clustering was used to establish groups de novo Then, probabilistic graphical models were used to study the differences between groups of patients at the biological process level. A molecular classification into two groups of patients was established, one group with increased expression of proteins related to adhesion, T lymphocytes and glycolysis; and the other group with increased expression of proteins related to translation and ribosomes. The functional analysis by the probabilistic graphical model showed that these two groups presented differences in metabolism, mitochondria, translation, splicing and adhesion processes. Additionally, these groups showed different frequencies of genetic variants in some genes, such as ATM, SLFN11 and DST Finally, genetic and proteomic characteristics of these groups suggested the use of some possible targeted therapies, such as PARP inhibitors or immunotherapy.

Genes and Variants Underlying Human Congenital Lactic Acidosis-From Genetics to Personalized Treatment.

Congenital lactic acidosis (CLA) is a rare condition in most instances due to a range of inborn errors of metabolism that result in defective mitochondrial function. Even though the implementation of next generation sequencing has been rapid, the diagnosis rate for this highly heterogeneous allelic condition remains low. The present work reports our group's experience of using a clinical/biochemical analysis system in conjunction with genetic findings that facilitates the taking of timely clinical decisions with minimum need for invasive procedures. The system's workflow combines different metabolomics datasets and phenotypic information with the results of clinical exome sequencing and/or RNA analysis. The system's use detected genetic variants in 64% of a cohort of 39 CLA-patients; these variants, 14 of which were novel, were found in 19 different nuclear and two mitochondrial genes. For patients with variants of unknown significance, the genetic analysis was combined with functional genetic and/or bioenergetics analyses in an attempt to detect pathogenicity. Our results warranted subsequent testing of antisense therapy to rescue the abnormal splicing in cultures of fibroblasts from a patient with a defective GFM1 gene. The discussed system facilitates the diagnosis of CLA by avoiding the need to use invasive techniques and increase our knowledge of the causes of this condition.

Draft Genome Sequences of Sporulation-Impaired Bacillus pumilus Strain NRS576 and Its Native Plasmid p576.

Bacillus pumilus spores can cause foodborne poisonings. B. pumilus strain NRS576 forms spores with a very reduced efficiency due to the presence of a plasmid, named p576. Here, we report the genome sequence of strain B. pumilus NRS576 and its plasmid p576.

Population structure of OXA-48-producing Klebsiella pneumoniae ST405 isolates during a hospital outbreak characterised by genomic typing.

OBJECTIVES: The aim of this study was to investigate the structure of a broad and sustained hospital outbreak of OXA-48-producing Klebsiella pneumoniae (KpO48) belonging to sequence type 405 (ST405). METHODS: Whole-genome sequencing and comparison of ten ST405 KpO48 isolates obtained from clinical samples in our hospital was performed. Using stringent criteria, 36 single nucleotide polymorphisms (SNPs) were detected (range 0-21 in pairwise comparisons), and allele-specific PCR was used to call the SNPs among a larger set of isolates. RESULTS: Several haplotypes were identified within the population. The haplotypes did not show a spatial structure, but a temporal evolution of sequential haplotype replacements was observed. CONCLUSIONS: The dispersed spatial distribution suggests a reservoir formed by a large pool of colonised patients, and the temporal replacement pattern suggests that the sustained outbreak was composed of several small outbreaks that appeared and rapidly dispersed to several units.

The transcriptional and mutational landscapes of lipid metabolism-related genes in colon cancer.

Metabolic alterations encountered in tumors are well recognized and considered as a hallmark of cancer. In addition to Warburg Effect, epidemiological and experimental studies support the crucial role of lipid metabolism in colorectal cancer (CRC). The overexpression of four lipid metabolism-related genes (ABCA1, ACSL1, AGPAT1 and SCD genes) has been proposed as prognostic marker of stage II CRC (ColoLipidGene signature). In order to explore in depth the transcriptomic and genomic scenarios of ABCA1, ACSL1, AGPAT1 and SCD genes, we performed a transcriptomic meta-analysis in more than one thousand CRC individuals. Additionally we analyzed their genomic coding sequence in 95 patients, to find variants that could orchestrate CRC prognosis. We found that genetic variant rs3071, located on SCD gene, defines a 9.77% of stage II CRC patients with high risk of death. Moreover, individuals with upregulation of ABCA1 and AGPAT1 expression have an increased risk of CRC recurrence, independently of tumor stage. ABCA1 emerges as one of the main contributors to signature's prognostic effect. Indeed, both high ABCA1 expression and presence of tumoral genetic variants located in ABCA1 coding region, seem to be associated with CRC risk of death. In addition the non-synonymous polymorphism rs2230808, located on ABCA1, is associated with gene expression. Patients carrying at least one copy of minor allele showed higher levels of ABCA1 expression than patients carrying homozygous major allele. This study broaden the prognostic value of ABCA1, ACSL1, AGPAT1 and SCD genes, independently of CRC tumor stage, leading to future precision medicine approaches and "omics"-guided therapies.